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Analytic approximants for the eigenvalues of the one-dimensional Schrodinger 
equation with potentials of the form V{x) = Ax"' + Bx^ are found using a multi- 
point quasi-rational approximation technique. This technique is based on the use 

b+2 

of the power series and asymptotic expansion of the eigenvalues in A = j4 a+2B, as 
well as the expansion at intermediate points. These expansions are found through 
a system of differential equations. The approximants found are valid and accurate 
for any value of A. As examples, the technique is applied to the quartic and sextic 
anharmonic oscillators. 



I. INTRODUCTION 



There are many techniques that have been developed with the purpose of solving the 
Schrodinger equation since it was first introduced more than eighty years ago. One might 
think that after so many years this should be a closed subject, yet this area of research is 
still being pursued nowadays by a number of physicists in the world. There are still many 
potentials of interest for which no exact solution is known and, moreover, new potentials 
that seem to model the behavior of physical systems, such as a set of different molecules 
interacting with each other or with an external field, are still being proposed from time 
to time, which prompts physicist to study them in more detail. A recent example of this 
can be found in Since one cannot find exact solutions for many of the most interesting 
potentials, one is left with one of two options: (1) use a numerical method or (2) try to find 
an approximate analytic solution. The last option is particularly appealing, since in many 
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situations it is possible to use precise approximate solutions in the same way as the exact 
ones. Approximate analytic solutions can be obtained by different methods for both, the 
energy eigenvalues and eigenf unctions of the Schrodinger equation, although more methods 
can be found in the literature for the eigenvalues than for the eigenfunctions. In either case, 
there are hundreds of publications devoted to this subject, and it would be impossible to 
make justice by citing them all. 

An analytic approximant should be a function of the parameters of the potential that 
comes very close to the values of the exact solutions (found numerically) when evaluated at 
any particular point in the parameter space. The usefulness of a particular method used to 
obtain these approximants will depend on the precision of the approximations as well as the 
simplicity of the analytic expressions. 

In this paper, a new method is proposed for finding analytic approximants to the energy 
eigenvalues of the one-dimensional Schrodinger equation with potentials of the form V{x) = 

. This kind of potentials are perhaps not the most interesting ones 
from a phenomenological point of view, but they are certainly of theoretical interest, since 
they include such potentials as the quartic and sextic anharmonic oscillators, and in any 
case, they are a good choice to test new methods before applying them to more interesting 
problems. 

The method we are proposing here is based on the multi-point quasi-rational approxi- 
mation technique, which has been successfull y appli ed before to similar kinds of problems 



involving differential equations j^, y, llO|, lUl, |12|, |13||- This technique consists in using the 



expansions of the function to be approximated around different values of the parameters in 
the differential equation where this function appears, in order to write an approximant in 
terms of rational functions in these parameters combined with auxiliary ones. The approx- 
imant will then have almost the same expansions around the different chosen values of the 
parameters. The auxiliary functions are usually needed in order to match the behavior of 
the quantity to be approximated when the parameters go to infinity, which normally cannot 
be done solely using rational functions as in a Pade's approximation. 

In our case, the differential equation we are interested in is, of course, the Schrodinger 
equation 

+ Ax" + Bx^]^ = Ei; . (1) 
ax^ J 

In this case, we have too parameters, A and B, which will be assumed to be both positive 
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to simplify the treatment, though the method can be used without this restriction. It is 
also assumed that a and h are positive integers, and b > a > 2. As it is very well known 
14| . one can make this equation to depend on only one parameter by making the changes, 

_ 1 _ 2 

X = A «+2x' and E' = A «+2£', which leads to 

(-^ + x- + Ax'^)v^ = £;V, (2) 

where A = A^°^B. From now on we will drop the primes and rename x' — > x and E' —>■ 
E. The energy eigenvalues E will depend on the parameter A. Our goal is to find an 
approximating function for E{X) for each energy level, using expansions around different 
values of A, including the power series (perturbative expansion around A = 0), and the 
asymptotic expansion (A oo). In sections HT] and UTT] a neat way to find these expansions 
will be shown using a system of coupled differential equations, which in the case of the power 
series provides an interesting alternative to standard perturbation methods, since here the 
perturbed eigenvalue can be found using only the corresponding unperturbed state, instead 
of using the whole unperturbed eigenvalue spectrum as in the usual quantum mechanics 
perturbation theory. In sections [IV] and |Vl the construction of the approximants will be 
shown, using the quartic and sextic anharmonic oscillators as examples. These approximants 
will be found for the ground state energy eigenvalues, as well as the first and second excited 
levels. The important point of our technique is that the same approximant will be valid and 
accurate for any value of A (including large and small values). 



II. POWER SERIES 

The expansion of the energy eigenvalues and eigenfunctions around A = can be written 

as 

E = E0 + E1X + EiX^ + ■ ■ ■ , (3) 
ijj = ^(, + ^^A + ?/'2A^ + -- - . (4) 

One would like to find the coefficients Eq, Ei, E2 . . .. This can be done introducing these 
expansions in equation ([2]), and demanding it to be satisfied at every order in A, which leads 
to the following system of differential equations 

Lipo = EqiPo , (5) 
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L^2 + 



(6) 
(7) 



n 



(8) 



where 



L 



.0, 



(9) 



It is important to note that, since A is arbitrary, many of the properties of the eigenfunction 
ip will be inherited by the functions ipoy V'l • • • • Fo^' example, given that for bound states, 
the function ip should fall off quickly for large values of x, so should the expansion functions 
ipk. Also, if the function ip has definite parity, then the functions ipk will all have the same 
parity as ip. 

The coefficients in the expansion of the energy can be found by solving numerically the 
differential equations one by one (using, for example, the shooting method). For instance, 
one could solve equation ([5]), obtaining then a numerical value for the coefficient Eq, together 
with a numerical solution for ipQ. For a bound state, one should find a solution for a range 
of values in x where the fall off to zero of ipo can be seen. Then one can make a very precise 
fit (a large polynomial in x) for ipQ within this range, and use this fit as an input to solve 
equation ([6]). This then leads to a numerical value for Ei, and a numerical solution for ipi, 
and in principle, the procedure can be repeated until one has as many coefficients E^ as one 
desires. 

Of course, the applicability of this method is limited by the numerical precision with 
which the differential equations are solved, and since the method is iterative, the numerical 
errors from the first n equations will be propagated to the solution of equation n + 1. For 
this reason, one would expect the precision of Ek to be lower for larger values of k. 

The numerical values of the coefficients E/^ can also be found directly using the solutions 
of the first k — 1 differential equations, without solving the k-th one. For example, if one 
has already obtained Eq and ipo{x), one can multiply both sides of equation ([6]) by ipoix), 
and integrating in x it is possible to show that 



El 



(10) 
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which coincides with the expression obtained using standard perturbation methods. The 
same procedure can be repeated for all the other equations, leading to 

This way of finding E^. is more precise, since the number of differential equation to be solved 
is — 1. In practice, the integrals are taken within a range of x where the functions ipkix) 
have already fallen to very small values. 

On the other hand, equation ([5]) can be solved exactly for a = 2, since then it would be 
the Schrodinger equation for a harmonic oscillator. It can be shown that in this case, all 
the other equations in the system can also be solved exactly. For example, for the ground 

2 

state Eq = 1 and ipo{^) oc exp(— ^). If we take 6 = 4 (quartic anharmonic oscillator) the 
next function, ipi{x), can be written as 

ijji (x) = {po + pix + p2x'^ + psx^ + PiX^) exp( — —) . (12) 

When this is introduced in equation ([6]), the function exp(— ^) disappears and a relation 
between two polynomials is left. Since this relation must be satisfied at each order in x, a 
system of equations in Ei and the p^'s is obtained, whose solution is 

3 13 
Pi = , p2 = -- , P3 = , Pa = --, Ei = -, (13) 

and it can be seen that po arbitrary, which means that just like for ipo{0), the initial condition 
■^1(0) = Po is arbitrary (this will be the case for all the other functions in the expansion). 
The same procedure can be repeated for 1^2, ips, etc., writing 

^n= (j^Pkxn expi-xy2) . (14) 



,fc=0 



We obtain 



^ 3 ^ 21 ^ 333 ^ 30885 , 

Eo = l, Ei = -, E2 = , E^ = , = . 15) 

' ^ 4 ' 16 ' 64 ' 1024 ^ ' 



This coincides with the results obtained by using the standard Rayleigh-Schrodinger per- 
turbation method, with the advantage that no information about the eigenstates of energy 
levels different from the one being considered is required in order to obtain the terms of 
higher order. The same can be done for other values of h. 
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Similar expansions can be found around any point other than A = 0. Let's call = A — a, 
then we can write 

+ + ax^ + Kx^ \i) = Ei) , (16) 



dx"^ 

and now we can expand around A^ = (i.e. around \ = a 



E = EX E'^K + E^Xl + ■ • • , (17) 
^ = ^o" + ^rA„ + 7/;°A2 +... . (18) 



The following set of equations is obtained 





LX 


= EX 




LX 


+ xVo 


= EX 


+ EX , 






= EX 


+ EX + EX 



where 



(19) 
(20) 
(21) 



LX + a^'C-i = E ^n-ki^^ ' (22) 

k=0 

L^ = -— + x^ + ax\ (23) 

Clearly, equation flTI?]) will not have exact solutions for any values of a and 6, so one is forced 
to find the coefficients numerically. Equations fllOj) and (ITT]) will still be valid, but with the 
changes iIju — > and E^ — > -E^ , i.e., 

r dxx\ij^f 

and 

The coefficient E'^ is actually the value of the fc-th derivative of the function E{X) evalu- 
ated at A = a. One might find these derivatives directly, evaluating E{\) nearby A = a, for 
example, 

= E(\ = <^)-E(X = a + ^) (26) 

However, this way of finding the coefficients becomes relatively difficult for higher derivatives, 
since then one needs to evaluate the function E{X) with increasing accuracy. The method 
proposed here can be viewed as an alternative, more accurate, easier and more efficient way 
to find these derivatives. 
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III. ASYMPTOTIC EXPANSION 

1 ~ 2+a ~ 2 

Doing the change of variables x = \ ^+''y, and defining A = A 2+i, and E = X the 



Schrodinger equation becomes 

rf2 



+ \y^ + y'']ij = E^ . (27) 



One can expand now E and t/; in a similar way as before, 

E = Ea + Ei\ + E2>? + ■■■ , (28) 

^ = ^^0 + ^13^ + W>^^ + ■ ■ ■ , (29) 

Introducing this in equation ( 1271) leads also to a system of differential equations 

L^o = ^0^0 , (30) 

LtAi + y>o = Eo^i + ^iV^o , (31) 

U2 + = ^0^2 + ^1^1 + ^2Vio , (32) 

n 

Lxpn + y^^n-l = ^ En-ki'k , (33) 



k=0 



where 



Rewriting the asymptotic expansion in terms of A instead of A, it is clear that the form 
of the expansion depends on the particular potential to be considered. In the case of the 
quartic anharmonic oscillator (a = 2 and 6 = 4), equation fl28l) leads to 

„ , 1 /'J / ^ E^ Eo E-i 

CO ^ oo -p^ _ oo 



= ^'"Eff + ^""E^ + lEl£^ (35) 

A:=0 k=0 k=0 

while in the case of the sextic anharmonic oscillator (a = 2 and 6 = 6), we obtain 



„ ^ 1 /A I E-\ Eo E-^ Ea 

A"*E# + ^-"*E%i- (36) 

fc=0 A:=0 
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In general, the expansions will have this structure, i.e., they can be divided in a few pieces, 
each one consisting in a series of negative integer powers of A, multiplied by a rational power 
of A. For this reason, the approximants that we will build must also be divided in a similar 
way, in order to match the behavior of each piece. This will be seen explicitly in the next 
two sections. 



IV. APPROXIMANTS FOR THE QUARTIC ANHARMONIC OSCILLATOR 



For the quartic anharmonic oscillator [15|, [la, ll7|, ll8|, ll9| , the approximants for the energy 



eigenvalues can be written in the following form 



where 



and 



N N N 



k=0 k=0 k=0 

N 



Q{X) = l + J2(lk^' ' (39) 

k=l 

that is, the approximant is constructed using rational functions multiplied by auxiliary 
ones, conveniently chosen in order to match the asymptotic behavior of the eigenvalues. 
Furthermore, since the power series is also going to be used, it should be possible to Taylor- 
expand these functions around positive values of A. It is for this last reason that the auxiliary 
functions are not chosen directly as the factors of A^^^, A^^^'^ and 1/A that appear multiplying 
each one of the three pieces that make up the asymptotic expansion. Instead, we do the 
change X —>■ 1 + fiX inside these roots, which, of course, still gives the right behavior for 
A — > oo. An arbitrary factor of /i has been included, which can be adjusted in order to 
improve the precision of the approximant. 

With this choice of auxiliary functions the degrees of the polynomials in the numerator 
must be the same as the ones in the denominator. In principle, this can be done indepen- 
dently for each one of the three pieces in equation fl37j) . i.e. Pa{X), Pb{X) and -Pc(A) could be 
chosen with different degrees, and in that case, different denominators matching the degree 
of each one of these polynomials would be needed. For simplicity, a denominator Q{X) com- 
mon to all three pieces has been chosen, and so all polynomials have the same degree. As it 
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will be understood later, any other choice would lead to a system of non-linear equations in 
the pfc's and g^'s, making the determination of the approximant unnecessarily complicated. 

The coefficients of the polynomials in the approximant are found using the power series, 
asymptotic expansion and the expansions around intermediate points (0 < a < oo), whose 
calculation was explained in the previous two sections. One is free to choose as many terms 
from each expansion as one desires, as long as the total number of terms from all expansions 
equals the total number of coefficients in the approximant. If the degree of the polynomials 
is A^, the total number of coefficients will be 4N + 3. In general, the approximants will have 
higher precision for higher N. 

The values of the first few terms in the power series (around A = 0) for the first three 
energy levels (labeled by n) of the quartic anharmonic oscillator are shown in table [H while 
the values of the first few terms in the asymptotic expansion are shown in table [Tll Notice 
that in accordance with what was discussed in section [Tll the values of the coefficients for 
the power series are exact. The values of the coefficients for the asymptotic expansion were 
obtained by solving equations (!30l) - (!33l) using the shooting method. For the ground state 
(n = 0) and the second excited level (n = 2), the eigenfunctions are even in x, and as 
mentioned before, so must be the functions ipk (and the same applies to ipk and so the 
initial conditions used in those cases were ipk{^) = 1 and ^'^(0) = 0. For the first excited 
level the eigenfunction is odd in x, so the conditions were V'fc(O) = and ip'j.iO) = 1. One 
might feel uneasy about the propagation of errors from one differential equation to the next, 
but it can be checked numerically that the accuracy of the energy eigenvalues for large values 
of A (or small values of A) improves as one includes higher terms in the expansion, which 
gives us confidence that the precision of the coefficients is acceptable. 

Let's choose a few intermediate points [i = 1,2,...), and let's take rij terms from 
the expansion around each one of these points. Let's also take uq terms from the power 
series (around A = 0) and Ua terms from the asymptotic expansion. It will be assumed that 
J2i ni + no + Ua = 4:N + 3. Using the power series at A = one can write 

g(A) V EkX"" = (1 + M)'/'^a(A) + (1 + fiX)-'/^Pb{X) + — — Pe(A) , (40) 

Taylor-expanding each side of this equation in A, and demanding it to be satisfied at every 
order up to A"°, one obtains a set of no linear equations in the coefficients of the approximant. 
Likewise, one can use the expansions at the intermediate points, and doing the change 
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Coefficients 


n = 


n = 1 


n = 2 


Eq 


1 


3 


5 


El 


3/4 


15/4 


39/4 


E2 


-21/16 


-165/16 


-615/16 


Es 


333/64 


3915/64 


20079/64 


E4 


-30885/1024 


-520485/1024 


-3576255/1024 


E5 


916731/4096 


21304485/4096 


191998593/4096 



TABLE I: Exact coefficients of tlie power series for tfie first three energy levels of the quartic 
anharmonic oscillator. 



Coefficients 


n = 


n = 1 


n = 2 


Eq 


1.060361944892 


3.7996728480480 


7.455697915983 


El 


0.362022935 


0.901605953 


1.244714261 


E2 


-0.034510565 


-0.057483095 


-0.046601602 


E3 


0.005195593 


0.005492673 


0.000958945 




-0.000831127 


-0.000513911 


-0.000831127 



TABLE II: Coefficients of the asymptotic expansion for the eigenvalues of the quartic anharmonic 
oscillator obtained solving the differential equations using the shooting method. 

A = Aa- + ttj, one can write 

N \ m N 



k=0 



fe=0 

N 

+(1 + /x(A„, + fti))-'/' Y ^^^^'-^ + "^)'' 

fc=0 

1 ^ 

+ (41) 



+ 



1 //(A„, + ai 



k=0 



If one demands this equation to be satisfied at every order in Aq. up to X^_, one obtains a 
set of rii linear equations in the coefficients. Finally, one can use the asymptotic expansion. 
For this we need to do the change X' — 1/X, and match the expansion with the approximant 
for each one of the three pieces in which it is divided. For example, since 

(42) 
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one can compare the term multiplying A^^^ in the right hand side of this equation with the 
term multiplying the same factor in the asymptotic expansion. Doing the same also for the 
other two pieces leads to 



iN-k 




N 



EskX"' = (A' + f^y/' Yl ' (43) 

k=0 k=0 

D N 

Y ^s.+iA'^'^ = (A' + ^^^""'^ ' (44) 



A:=0 k=0 

oo N 

fc=0 ^ k=0 

Here the number of terms taken in each expansion is determined by n^, that is, one would 
not allow any Ek with k > Ua itl the sums. In this way, one gets a set of Ua linear equations 
for the coefficients of the approximant. 

In table IIVI the values of the coefficients of the approximants are shown for the first 
three energy levels, using polynomials of degree three. There are fifteen coefficients in 
each approximant, and they were obtained using the first five terms of the power series 
(around A = 0), the first five terms of the asymptotic expansion, and the first term of the 
series around A = 0.5, A = l, A = 2, A = 5 and A = 20 (which are shown for the three 
energy levels in table HTB . This means that we are only using the exact energy eigenvalue 
around these intermediate points, and forcing the approximant built with the power series 
and asymptotic expansion to furthermore coincide with these "exact" eigenvalues at these 
points. This not only brings the relative error of the approximant at these points down to 
zero (they become nodes of the relative error as a function of A), but also helps to decrease 
the error in between these points. The relative error is defined using as target the eigenvalues 
obtained numerically through the shooting method, i.e., the relative error is given by 

I Eapp — -^shooting I ^^g^ 
-^'shooting 

The highest relative error with these approximants was obtained for small values of A. 
Specifically, the maximum relative error was obtained around A ~ 0.2. In the case of the 
ground state, the highest relative error was 

I Eg^pp -£'shooting I 



-E^shooting 



= 1.05 X 10"'' . (47) 

A=0.17 

The relative error decreases rapidly for smaller values of A, and of course, it also decreases 
when A increases until it finds the next node at A = 0.5. After that, the relative error never 
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CoefRcients 


n = 


n = 1 


n = 2 




= 1/2) 


1.241854043136 


4.051932338617 


7.396900686938 


Eo{X 


= 1) 


1.392351580103 


4.64881282723 


8.655049982254 


Eo{X 


= 2) 


1.607541348124 


5.475784646286 


10.358583364736 


Eo{X 


= 5) 


2.018340657447 


7.013479298703 


13.467730394819 


Eo{X 


= 20) 


3.009944947791 


10.643215959124 


20.694110927154 



TABLE III: First coefficient (energy eigenvalues) of the series at different intermediate points for 
the first three energy levels with V{x) = + Xx^. These values were obtained using the shooting 
method. 

becomes higher than 2 x 10~^. For the first and second excited level, the maximum error 
around A ~ 0.2, was about 8 x 10"'' and 2.4 x 10"^, respectively, and after the node at 
A = 0.5 this error is never higher than 4 x 10"*^. In fact, the relative error decreases quite 
rapidly in the case of the first and second excited levels for large values of A, although it 
does so more slowly in the case of the ground state. 

In all these approximants, we chose n = 2. This parameter is arbitrary except for one 
restriction: The approximants should not have any defects, that is, there should not be 
any poles in the approximant (positive roots of (5(A)) with the corresponding nearby zeros. 
Notice in table [IV]that with this choice of /i all of the coefficients of Q{X) are positive, which 
will, of course, guarantee that it has no roots for A > 0. Other choices of fi may lead to 
mixed negative and positive coefficients in Q{X), which will in general lead to positive real 
roots in this polynomial. Other than that, there is no other restriction in fi. Among all of 
the values of fi that allow to keep the approximant free of defects, one is free to choose the 
one that minimizes the relative errors. 

Other than improving the numerical method used to obtain the coefficients of the ex- 
pansions, there are several ways in which the maximum relative error of the approximants 
can be decreased for all energy levels. The easiest one is to move one of the nodes. For 
example, one may choose the approximant to have a node at A = 0.2 instead of A = 0.5. 
If this is done, it can be seen that the maximum relative error is reduced by about a half 
for all energy levels. Another possibility is to use the derivatives of E{X) at some of the 
intermediate points. Finally, one may try an approximant of higher degree, allowing it to 
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Coefficients 


n = 


n = 1 


n = 2 


ao 


-235.587774594 


3.26113271857 


-46.3903727540 


ai 


129.192528081 


45.7861842084 


118.622015136 


a2 


819.219808968 


347.592601172 


906.778841942 


as 


4083.20247083 


1023.55148495 


3353.40199807 


bo 


49.9309955808 


4.80222464913 


8.35806073416 


bi 


374.551951382 


43.2593054339 


113.555378592 


62 


1181.63222463 


259.439145729 


536.540936096 


bs 


2212.93937770 


385.537761596 


888.696857827 


Co 


186.656779013 


-5.06335736770 


43.0323120198 


Cl 


208.866219507 


-20.7666223608 


48.6886778830 


C2 


-423.418335743 


-35.6233569984 


-59.8483255497 


C3 


-334.866518168 


-39.0190739652 


-52.8167275564 


91 


148.201294158 


24.5427291322 


29.7104983824 


92 


1782.00574019 


171.823102857 


247.681714807 


93 


4851.65727491 


339.396077796 


566.683604101 



TABLE IV: Coefficients for the approximants of the ffi'st three energy eigenvalues of V{x) = 
+ Ax^, using polynomials of degree 3. 

coincide with the values of E{X) and its derivatives at more points. This last possibility is 
studied in the next example. 



V. APPROXIMANTS FOR THE SEXTIC ANHARMONIC OSCILLATOR 



For the sextic anharmonic oscillator 20|, l2l|], the asymptotic expansion consists of only 



two pieces, so the approximant can be written as 



£;app(A) = (1 + ^^>^Y"^^ + (1 + M)" 



1/4 



n(A) 



(48) 



g(A) ^ " ' Q{\) ' 

where -Pa(A), -Pfe(A) and Q{\) are given in equations (!38l) and (!39|) . The corresponding 
coefficients a/c, hk and are found in a similar way as we did for the quartic anharmonic 
oscillator. For approximants of degree iV, we will have 3A^ + 2 coefficients, so taking tiq 
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Coefficients 


n = 


n = 1 


n = 2 


Eq 


1 


3 


5 


El 


15/8 


105/8 


375/8 


E2 


-3495/128 


-47145/128 


-295095/128 


E3 


1239675/1024 


27817125/1024 


276931275/1024 


Ea 


-3342323355/32768 


-110913018405/32768 


-1626954534555/32768 



TABLE V: Exact coefficients of the power series for the first three energy levels of V{x) = + Xx^. 

terms from the power series, rij terms from the series at the i-th intermediate point, and 
terms from the asymptotic expansion, we should have + no + = 3A^ + 2, and the 

coefficients a^, and can be determined using equations derived by matching powers in 
A in the equation 

no 



Q{\) J2 EkX' = (1 + M)i/^P«(A) + (1 + fiX)-'/'P,{X) , (49) 

fc=0 

matching powers in Aq,- using 

N \ rii N 

l + Y.qdXa.+a,r]Y.ET^^. = (l+MA., + «.))^/^5^a.(A„,+a.)^ 

k=0 

N 

+(1 + /i(A„, + a,))-i/^ Yl ^'^(^"^ + «^)' ' 



k=l 



k=0 



k=0 



and matching powers in A' in 



N 



V k=l / k=0 k=0 

N \ 00 N 



iN-k 



k=l 



k=0 



k=0 



(50) 

(51) 
(52) 



The first few coefficients of the power series at A = and the asymptotic expansion, 
obtained using the systems of differential equations described in sections Hi] and IIIII are 
shown in tables |V] and IVTI respectively. As expected, the coefficients of the power series are 
exact, and the coefficients of the asymptotic expansion are obtained numerically. 

The degree of the polynomials in the approximant was first chosen to be = 5. The co- 
efficients of the approximants were calculated choosing different intermediate points (nodes) 
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Coefficients 


n = 


n = 1 


n = 2 


Eq 


1.144802430723 


4.338598612643 


9.073084583078 


El 


0.307920324 


0.718220191 


0.904435602 


E2 


-0.018541674 


-0.024395762 


-0.010249211 


E3 


0.001559745 


0.00099946565 


-0.000749318 


Ei 


-0.000123969 


-0.000026329 


0.000107946 



TABLE VI: Coefficients of the asymptotic expansion obtained solving the differential equations 
using the shooting method for V{x) = + . 

for the first three energy levels, together with the first four terms of the power series and 
the first five terms of the asymptotic expansion. For the ground state, the intermediate 
points were A = 0.1, A = 0.2, A = 0.5, A = l, A = 2, A = 5 and A = 10, and only the 
energy eigenvalues at those point were used, i.e., we didn't pick any of the derivatives at 
these points. Furthermore, we chose /i = 1/2. For the first and second excited level, we 
used El, El, El, E^^, E^^, E^-^ and with fi = 0.95. Notice that in these cases, 

the derivative at A = 0.1, i.e., Ef'^, has also been used. With these choices, the values of 
the coefficients in the approximants are given in table IVIIi With these approximants, the 
highest relative errors were obtained around small values of A. Specifically, at A = 0.014 the 
relative error for the ground state eigenvalue is 2.55 x 10~^. For even smaller values of A, 
this error decreases rapidly. For 1 < A < 5 that maximum error is 1.3 x 10~^, and for A > 5 
the maximum error was 7 x 10~^, decreasing rapidly for large values of A. 

As it can be seen, in both the quartic and sextic anharmonic oscillators, the region of A 
where it is more difficult to achieve high accuracy is for A < 0.5. We tried also approximants 
of degree 6 for the sextic anharmonic oscillator, and it was possible to reduce the relative 
error in this region by a factor of 1/2. The coefficients of the approximants can be found 
in table IVIIII For n = 0, 2, we used /i = 1/2, while for n = 1, we used /i = 1. For n = 
the approximant was built using the first four terms of the power series around A = and 
the first five terms of the asymptotic series, together with E^, Eq, Eq, Eq^, E^^, E^"^ , E\^'^ , 
-Eg'^, -Eq'^, -^i'^ and Eq'^^. For n = 1,2, the same terms were taken, except that the fourth 

1 /2 

term of the power series around A = was replaced by E2 (i.e. the second derivative at 
A = 1/2). The improvement in the accuracy of the approximants for A < 0.5 is related to 
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Coefficients 


n = 


n = 1 


n = 2 


Oo 


-228343.425175410234 


-1455854.05235538211 


-1636184.16769173318 


Ol 


-51504.1866802877753 


-523067.617460085793 


-319382.170585718291 


02 


-3246.30836836582236 


-82992.9062424631909 


-76246.2680853710652 


03 


17641.2930769775453 


564652.836765077742 


1736863.85800628006 


04 


40455.7925641666977 


2737870.84947363967 


8425191.38724236837 


05 


25845.5374941939783 


3765927.01323457233 


10377310.4646700931 


^0 


228344.425175410234 


1455857.05235538211 


1636189.16769173318 


^'i 


108718.883825186570 


1215348.76068812833 


1098149.56754718475 


b2 


12471.9577982789445 


211783.949326608282 


156879.990271492461 


bs 


8971.90069757654071 


247389.200455627253 


461008.702437356268 


64 


12714.5141045106691 


765825.922719107915 


1349548.74023090975 




4915.62090727820199 


607633.680665566723 


1008252.48843727777 


9i 


126.840851046236208 


245.543619745344247 


306.370961578640294 


Q2 


3258.74684448027568 


13846.7357289482276 


20215.7966313260916 


Q3 


21339.1202042371419 


208223.456783358953 


314306.763606766761 


Q4 


39515.8524539180998 


853339.703125454811 


1215186.64688146943 


Q5 


18984.4284445198520 


856945.745673868394 


1129173.69341497105 



TABLE VII: Coefficients for the approximants of the first three energy eigenvalues of V{x) 
+ Ax^, using polynomials of degree 5. 



having used more terms in the expansions around points in this region, including first and 
second derivatives. 

The fact that the approximants become less accurate for small values of A may be related 
to the analytic properties of the exact function E{\) [l^, although this is far from clear. As 
it is well known, the perturbative expansion (i.e., the power series around A = 0) is divergent 
for any A 7^ 0. However, in previous works where quasi-rational approximants have been 
used, it has been found that the main factor determining the accuracy of an approximant is 
;he accuracy of the coefficients in the expansions used to build it. For example, in references 



22] and 23j series with radius of convergence equal to zero were used, yet the approximants 



were very accurate because the coefficients of these series were determined with very high 
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precision. 
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Q3 


157632.217425587120 


133763.388618245199 


150748.225437877810 


Q4 


775120.602962958890 


586955.795810181245 


587462.218227788458 


Q5 


1249527.64792219401 


727254.753026886137 


723876.015243016534 


96 


564727.576025957156 


102119.617954584846 


248600.057576103110 



TABLE VIII: Coefficients for tfie approximants of the first tfiree energy eigenvalues of V{x) = 
Xx^, using polynomials of degree 6. 

VI. CONCLUSIONS 

In this paper, it has been shown that accurate approximants for the energy eigenvalues of 
potentials of the form V{x) — Ax"' + Bx^ can be found using a multi-point quasi-rational ap- 
proximation technique. The approximants are constructed using rational functions, together 
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with auxiliary functions introduced to be able to reproduce the behavior of the eigenval- 
ues for large A. The coefficients of the rational functions are found using the power series 
of the eigenvalues, not only at A = 0, but also for arbitrary finite values of A, together 
with the asymptotic expansion. These expansions are found using a system of differential 
equations, which in the case of the power series at A = 0, represents an alternative way to 
find the perturbative expansion. As examples, approximants for the lowest energy levels of 
the quartic and sextic anharmonic oscillators were obtained. The approximants were fairly 
simple, since the degree of the polynomials used was not too high. In particular, for the 
quartic anharmonic oscillator, it was shown that it is possible to obtain approximants with 
polynomials of degree 3 for which the relative error is not higher than ~ 3 x 10~^. In the 
case of the sextic anharmonic, polynomials of degree 5 and 6 were tried, and it was shown 
that in the second case it was possible to improve the accuracy for A < 0.5 while maintaining 
the accuracy for larger values of A 

In this technique, one has a lot of freedom choosing the intermediate points, as well 
as how many terms from each one of the series to take. This gives a lot of possibilities 
to try to reduce relative errors, and it would be interesting to study if it is possible to 
do this systematically. It would also be interesting to try to find methods that allow to 
improve the accuracy of the coefficients in the expansions. For example, if better numerical 
solutions for the expansion functions can be found, this should lead to better coefficients 
and therefore to better approximants. As it was mentioned before, experience with quasi- 
rational approximants has shown that the accuracy of the coefficients in the expansions is 
the main factor infiuencing the accuracy of multi-point quasi-rational approximations. For 
simplicity, here we have limited ourselves to the use of the shooting method to solve the 
systems of differential equations. On the other hand, the first equation in each one of the 
systems, i.e., equations ([5]), (fT9|l and (l30!l are just regular Schrodinger equations, and many 
methods have been developed that allow to obtain the corresponding eigenvalues with very 

nnnnnnnnnnn 

high precision [5|, 117|, 118|, 119|, Um, 1251, 126|, |27|, |28!, 1291, 130|. In fact, we could have used any of 



these methods to obtain the coefficients Eq and Eo, but it is not clear how to extend these 
methods to solve the remaning equations in the systems. In the future, we plan to study 
these issues in more detail, and apply this technique to other potentials of interest, both in 
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physics as well as chemistry. 
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